ON 



NORDITA-96/18 P 
March 1996 
hep-ph/9603266 



Color-Octet J/i/j Production at Low p± 

NO ■ Wai-Keung Tang9 

on ; 

Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



M. VanttinenUl 
NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 
(February 1, 2008) 

> ' 

no 

NO 
(N 
m . 

O . Abstract 

NO 

On 

& ' We study contributions from color-octet quarkonium formation mechanisms 

O j' to J/ij) hadroproduction at low p±. We include transitions of color-octet cc 

^ ' states into "direct" J/ip and into xi,2 which decay radiatively into a J/tp. 

^ ■ Together with earlier work, this calculation constitutes a complete analysis of 

p^-integrated J/ip production at leading twist. We find that the leading-twist 
contribution is not sufficient to reproduce the observed production rates and 
polarization of the J/ip and Xi,2- Hence there must exist other important 
quarkonium production mechanisms at low p±. 
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I. INTRODUCTION 



The production and decays of heavy quarkonia have been widely studied using perturba- 
tive QCD M. Due to the large masses of the c and b quarks, the production or annihilation 
of heavy- quark-ant iquark pairs takes place at a much shorter timescale than the forma- 
tion of the bound state. This makes it possible to factorize the transition amplitudes. In 
hadroproduction reactions, also the initial-state hadronic structure is usually expected to 
be factorizable from the QQ production dynamics in terms of single parton distributions; in 
other words, quarkonium production is usually taken to be a leading-twist process. 

In the simplest approach, one takes the hadroproduction cross section to be a fixed 
fraction of the integrated QQ cross section below the open heavy flavour threshold. This 
is called the semilocal duality model or color evaporation model. This model reproduces 
successfully the dependence of the cross sections on the center-of-mass energy of the colli- 
sion and on the longitudinal momentum of the quarkonium 0. However, many observables 
such as the relative production rates of different charmonium states, transverse momentum 
distributions, and quarkonium polarization are not predicted. In a more sophisticated ap- 
proach, a color-singlet QQ pair with the appropriate J PC quantum numbers is produced in 
the hard process. The quarkonium production amplitude is then written as the convolution 
of a perturbative amplitude for the production of a QQ 2S+1 L^P pair (the superscript "(1)" 
stands for a color-singlet configuration) and a nonrelativistic wave function. This model, in 
which only the leading, color-singlet Fock component of the quarkonium wave function is 
taken into account, is called the color singlet model or charmonium model. 

In the color singlet model, the polarization of final-state quarkonium is determined by 
the perturbative dynamics of QQ production and by the angular momentum projection in 
the wave function. Since polarization is relatively insensitive to higher-order corrections, it 
provides a very good probe of the basic QQ production mechanisms. 

Experimentally, J/ip polarization has been measured in fixed-target tc~N 0-§J and pN 
reactions. The parameter a in the angular distribution 1 + a cos 2 9 of J/ip decay dileptons 



2 



in the Gottfried- Jackson frame has been measured to be a = 0.028 ± 0.004 in ir~N reactions 
and a = —0.115 ± 0.061 in pN reactions for longitudinal momentum fraction xp > at 125 
GeV ||. This corresponds to unpolarized production. The theoretical interpretation of 
these results is complicated by the fact that the "direct" and Xc-decay components of J/ip 
production have not been resolved in the polarization analysis. However, the polarization 
of the ip'(2S) has been measured in tt~N reactions ||. The if}' is also produced unpolarized, 
with a = 0.02 ± 0.14 for xp > 0.25 at 253 GeV. Apart from an overall normalization factor, 
the direct J/i[) cross section is expected to be similar to the total ip' cross section, because 
no significant contribution to ip' production from the decays of higher-mass states has been 
observed. We therefore assume that the direct J/ip component is unpolarized; the x-decay 
component should then also be unpolarized. 

In Ref . |7j , we calculated the leading-twist contribution to the p_i_-integrated cross section 
<j(-k~ N — > J/ip(\) + X), where A is the helicity of the J/ip, within the color singlet model. 
We included contributions from the direct mechanism gg —>■ J ftp + g and from the radiative 
decays of xj states produced in the reactions gg — > X2, 99 —> Xi9i 19 ~^ Xili an d qq — > Xi9- 
The polarization analysis of the decay contribution is made possible by the electric dipole 
nature of the decay. In ip' production, only the direct mechanism needs to be taken into 
account. 

The leading-twist, color-singlet contributions to J/ip and ip' production turn out to be 
dominantly transversely polarized, i.e. a(X = ±1) > er(A = 0), with a ~ 0.5 for the total 
J/ijj and a ~ 0.3 for the ip' and direct J/ip produced by pions. This is in contrast to the 
experimental observation of unpolarized production. The model also fails to reproduce the 
observed relative production rates of the J/ip, Xi an d Xi states. In Ref. 0, we inter- 
preted these discrepancies as evidence for important higher-twist mechanisms of charmonium 
production. 

However, one could argue that the discrepancies are due to neglecting the higher Fock 
components in the quarkonium wavefunction. These are systematically included in nonrel- 
ativistic QCD (NRQCD), an effective field theory which has been formulated during the 
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recent years p|,|TO||. NRQCD provides an expansion of quarkonium cross sections and decay 
widths in terms of the relative velocity v of the heavy quark and antiquark in the bound 
state. The nonperturbative physics is factorized into an infinite number of matrix elements 
that scale in a well-defined way in rriQ and v. 

Color-octet mechanisms have recently been suggested jTTHTBI] as an explanation of the 
large quarkonium production rates observed at the Tevatron pp collider |l4lJl5fl . Some of 
these rates exceed color-singlet-model predictions by more than an order of magnitude. 

There exist color-octet mechanisms of low-p_i_ quarkonium production which are of higher 
order in v 2 but lower order in a s compared to the leading color-singlet mechanisms. They 
might help to reproduce the fixed-target data; the purpose of this paper is to find out 
whether they do. If not, then the fixed-target discrepancies must be due to a breakdown of 
the leading-twist approximation of QQ production. 



In a recent work | 16| , we already analyzed the color-octet contribution to if)' hadroproduc- 
tion at low p±. We found that the leading color-octet contribution is dominantly transversely 
polarized. Hence, even with the leading color-octet components included, the unpolarized 
if)' production data || cannot be reproduced. This suggests that there are important higher- 
twist QQ production mechanisms. The ip' analysis also applies to the direct component of 
J I if) production. In the present paper, we complete our analysis by calculating the contribu- 
tion to J I if) production from the radiative decays of \i an d X2 produced through color-octet 
intermediate states. It will turn out that this component, too, is dominantly transversely 
polarized. Furthermore, determinations of color-octet matrix elements from the analysis of 
other reactions imply that the magnitude of this component is small compared to the data. 
We shall conclude that in J/ if) as well as in if)' production, new mechanisms are likely to be 
important at low 
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II. SUMMARY OF COLOR-SINGLET J/tp PRODUCTION CROSS SECTIONS 



At leading order in ot s and at leading twist, the color-singlet QQ production subprocesses 

are 

gg 'So, 3 P ,2, (2.1) 
gg^^ + g, 3 Pj + g, (2.2) 
gq^ 3 Pj + q, (2.3) 
qq -> 3 Pj + g, (2.4) 

i.e. they are obtained from the leading-order annihilation subprocesses such as J/ip — > ggg 
by crossing. Note that the process gg <-> 3 Pi, with the gluons on mass shell, is forbidden by 



Yang's theorem ||17|| . Within the color singlet model, the QQ production mechanism could 



be different from (|2.1| - l2.4|) ; it could e.g. be a higher-twist process, where many partons from 
the same initial-state hadron participate in the hard scattering. 

The radiative decays Xi,2 ~> J/ip + 7 are known experimentally to be a major source 
of J/ip production. They contribute 30-40% of the total cross section in both pion- and 
proton-induced reactions ||. The production of the ip', on the other hand, is expected to 
be dominantly due to the direct subprocess (|2.2j). 

The polarization of a 3 S\ state like the J/ip is reflected in the polar-angle distribution 
of its decay dileptons in their rest frame. The parameter a in the angular distribution 
1 + a cos 2 6 is related to the polarized J/ ip production cross section 

o-(A) = fftot -r- , 2.5 

6 + a 

where A is the helicity of the J/ip, in the following way: 

_ da{\ = 1) -2dcx(X = 0) + da{\ = -1) a 
a ~ da{\ = 1) + 2da{\ = 0) + da{\ = -1) ~ ~2 + a 

Because of rotational and parity invariance, da(X = 1) = da(X = —1). We work in the 
Gottfried- Jackson frame, which is defined as the particular quarkonium rest frame where 



the beam momentum lies along the z axis and the target momentum lies in the xz plane, 
with p x (target) < 0. At vanishing transverse momentum of the quarkonium, the Gottfried- 
Jackson frame is obtained from the laboratory frame by a simple boost, and the choice of 
the x axis is immaterial. 

We analyzed the polarized cross sections of J/ip production through the direct process 
and xj decays in Ref. [0, where we found that the color singlet model is insufficient to 
explain the existing fixed-target data. In the color singlet model, the xi an d direct J /if} 
cross sections are predicted to be significantly lower than measured relative to the xi cross 
section, which is within a factor K = 2-3 of the experimental value. In contrast with 
the observed unpolarized production, the predicted total J /if) cross section is dominantly 
transversely polarized even if the various contributions are renormalized according to the 
data. The quantitative results of the calculation of Ref. are listed in Table | together 
with experimental data and color-octet predictions. 



III. COLOR-OCTET J/ip PRODUCTION CROSS SECTIONS 

As one calculates quarkonium production cross sections within the NRQCD factorization 
scheme, one is using two expansions: the perturbative expansion of the short-distance QQ 
production amplitude and the velocity expansion of the long-distance quarkonium formation 
amplitude. General rules for finding out the power dependence of the NRQCD matrix 
elements on rriQ and v can be found in Refs. 0|lOj. At leading order in perturbation theory, 
i.e. O(o%), and up to next-to- leading order in the velocity expansion, the subprocesses for 
leading- twist J /if) production through color-octet intermediate states are 

J/if> + gg, (3.1) 

J/4 + 9, (3.2) 

J/iP + g, (3.3) 

XJ + 9 -> J/^ + 7 + 9- (3-4) 
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These are illustrated by the Feynman diagrams of Fig. [I], where the blob represents a nonper- 
turbative transition. Their cross sections are proportional to the NRQCD matrix elements 

(3.5) 
(3.6) 
(3.7) 
(3.8) 

respectively. Hence the direct J/ijj production cross sections are proportional to a 2 v 7 , and 
the xj cross section is proportional to a 2 v 5 . They are to be compared with the leading 
color-singlet cross sections 

o{gg -> J/ip + g) ~ «^ 3 , 

a(ij -> xi + *:) ~ "> 5 . (3.9) 



The amplitude 



qq 



cc 



3p(8) 



J 



■J/ip + g 



(3.10) 



is of higher order in v 2 than the amplitudes of the processes (|3.1| -|3~4l), because the lowest- 
order non-perturbative transition (single chromoelectric dipole) is forbidden by charge con- 
jugation. Furthermore, the amplitude 



gg 



cc 



3g(8) 



Xj 



(3.11) 



is of higher order in a s because the amplitude A (gg — > cc 3 S'i 8 ' ) J vanishes in the leading 
order. Charge conjugation or Yang's theorem would not require this amplitude to vanish 
because the cc pair is not in a color-singlet state. 

The contribution to a(hN — > J/i[)(X) + X) from the direct production subprocesses ( J3. 1| — 
follows immediately from the ip' production analysis of Ref . jTjJ : 
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VoctetihN -> direct J/V(A) + X) = O 1 (0|O 8 J/V '( 3 Pi)|0>(3 - 25 



AO J 

J/i>/ 1 , 



+ o 2 (o|o 8 ^(%)|o)(i-M 

+ O 3 (0|O 8 7/ ^( 3 5 1 )|0)(l-M- (3-12) 



The coefficients are 



Oi = "g^r $ gg /^(M 2 / S , axjp), (3.13) 

° 2= 24Mt $ ^ Jv(M2/s ' /iF) ' (3 - M) 
03 = §5 E [«Wwv(M7s, //f) + ^ q/hN (M 2 /s^ F )} , (3.15) 



where 



$ij/hN{T- 



fjip) = J dx 1 dx 2 fi/h(%i, Vf) fj/N(x 2 , Hf) (3.16) 

is a parton flux factor evaluated at the leading-twist factorization scale /ip (since the final- 
state gluons are taken to be soft, the kinematics is essentially that of a 2 — > 1 subprocess, 
with s = M 2 = 4m 2 ). For ease of reference, we plot in Fig. ^ the gg and qq flux factors 
in 7r~-proton, proton-proton and antiproton-proton collisions using the GRV-LO parton 
distributions [18]. In accordance with the existing experiments, the integral is over the 
region x F = X\ — x 2 > 0. 

The polarization of the directly produced J / ip has not been measured separately from the 
polarization of those from x decays. However, the fact that the ip' are produced unpolarized 
H lets us expect that the direct J / tp component is also unpolarized. Then, as in the if)' case, 
only about half of the observed direct production can be due to the strongly transversely 
polarized contribution (|3.12|) . We discussed the derivation of a quantitative bound on a 
linear combination of NRQCD matrix elements in Ref. ||16|| . 



To evaluate the contribution from the process ( |3.4j ), we make use of the fact that both 
the cc 3 S[ 8 ^ — > xj + 9 transition and the radiative decay of the xj are electric dipole 
transitions. The necessary formulas are given in the appendix. From the qq — > cc scattering 
amplitude, 
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The subprocess cross sections then become 
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(3.17) 



(3.18) 



(3.19) 



(3.20) 



Using the NRQCD result (0|C»f ( 3 Si)|0) = (5/3) (O|0f ( 3 Si)|0), valid up to corrections of 
relative order v 2 , we arrive at 



= E 

3 
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J/V> + 7) 



47-215 



AO 



120 



(3.21) 



Our expressions for the polarized cross sections, eqs. (|3.12| , |3.21| ), agree with the unpolarized 
cross sections given in Refs. [lT3] , p~9| . 



Each of the two components in eq. ( |3.21| ) alone would correspond to a = 0.2 and a = 0.29 
in the decay angular distribution of the J/ip, respectively. Hence this contribution, similarly 
to the color-singlet contribution M and the direct color-octet contribution of eq. (|3.12j ), is 
dominantly transversely polarized. Because J/ip production in irN reactions is observed 
to be unpolarized, one is again led to conclude that there are important cc production 
mechanisms beyond leading twist. 

The value of the (0|Og 1 ( 3 jS'i)|0) matrix element has been obtained from the analysis of 
CDF data at large p± [|T2"| , |i"3l and of B decays into P-wave charmonium [[TI|,[2(J. These 



values actually imply that the contribution ( |3.21|) is rather small compared to the observed 
J/ip cross section. Using (0|Of ( 3 Si)|0) = (9.8 ± 1.3) ■ 1(T 3 (GeV) 3 @, a s = 0.26 and 
M = 3.5 GeV, we obtain 



J2 a octct(hN -> xi,2 + 9 -> J/^(A) + 7 + fl0 

A 

= 4.5 nb 



/ M 2 \ / M 2 N 

®qq/hN , A*F + ®qq/hN , A*F 



1.3 nb (7r beam) 
0.6 nb (p beam) 



(3.22) 



at £^i a b(7r) = 300 GeV. These numbers are more than an order of magnitude smaller than 
the experimental cross sections of 72 nb, 67 nb and 45 nb (with errors of about 25%) for 7r + , 
7r~ and p beams, respectively ||. Hence the color-octet mechanisms cannot reproduce the 
fixed-target data. On the other hand, our analysis does not set any constraints that would 
contradict the color-octet description of other reactions. 



IV. SUMMARY 

In this paper, we have considered the production of J/ip charmonium in fixed-target 
reactions within the factorization scheme of nonrelativistic QCD (NRQCD) ||10|| . We have 
calculated the contribution from the radiative decays, Xi,2 — > J/ip + J, of Xj charmonia 
produced through intermediate color-octet cc states at leading twist. The "direct" color- 
octet component of J/ip production is given by our analysis of color-octet ip' production 
16| . Together with our earlier evaluation of the color-singlet component [f5]], the present 



calculation constitutes a complete analysis of leading-twist J/ip production at low p±. We 
have included contributions from nonperturbative transitions between intermediate cc states 
and charmonium up to relative order t> 4 , where v is the relative velocity of the charm quark 
and antiquark in the bound state. The cc production amplitudes have in each case been 
evaluated at the lowest order in a s allowed by the quantum numbers of the intermediate cc 
state. 

The results on J/ ip production obtained in this paper and Refs. |7|JT6| have been collected 
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in Table |I[ In the evaluation of color-octet contributions, we have used the NRQCD matrix 



elements determined in Refs. 13.21 



Our motivation has been to test both the NRQCD picture of the long-distance quarko- 
nium formation process and the leading-twist approximation of short-distance QQ produc- 
tion. Predictions of the final-state charmonium polarization provide a very good test of the 
models involved. They are relatively insensitive to higher-order corrections in perturbation 
theory. Polarization analysis of the long-distance process is made possible by the fact that 
the emission or absorption of soft gluons does not change the heavy quark spins. 

We have shown that the leading-twist contribution to J/tp production is dominantly 
transversely polarized, i.e. a > in the angular distribution, 1 + a cos 2 9, of the decay 
J/ip — > in the Gottfried- Jackson frame. Actually, all but one of the individual com- 
ponents of the theoretical cross section are dominantly transversely polarized. The only 
exception is the small contribution from the radiative decay of the xi produced by color- 
singlet mechanisms, which gives a ~ —0.15. Experimentally, on the other hand, it has been 
observed that S"-wave charmonia are produced unpolarized. 

Furthermore, the values of NRQCD matrix elements determined from other reactions 
imply that the normalization of color-octet contributions is small compared to the observed 
charmonium cross sections. 

Assuming that the NRQCD factorization scheme provides a complete description of 
quarkonium formation, the reason for the discrepancies must lie in the choice of cc production 
mechanisms. Hence there should exist important higher-twist mechanisms of cc production 
at small p±. 

Since the mass of the b quark is significantly larger than the c quark mass, all the approx- 
imations involved in the calculation - perturbation theory, the velocity expansion, and the 
leading-twist approximation - are expected to work better for bottomonium. Unfortunately, 
the existing bottomonium production data is insufficient to test our predictions. 
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APPENDIX A: DIPOLE TRANSITIONS 



1. Electric dipole transitions 

The polarized cross section of J/ip production via the radiative decay of a xj charmonium 
state is 

= 1 J dLips( J/rl>, 7, Xf£ \A (ij ^ XJ + X^ J/^(\) + 7 + X) | 2 

= ^ / dLips( X j, X) dLips( J/^, 7 ) EH 2 > ( A1 ) 



where s = (pi + p.,) 2 is the invariant mass of the initial state and A is the helicity of the J ftp. 
In the limit where the momenta of the photon and final-state light hadrons are neglected, 
it equals the z component of the spin of the J/ip. The transition amplitude is 

A (ij -^ XJ + X^ J/i>(\) + 7 (/i) + X) 

= Y,A (ij - X j(Jz) + X) v( Xj ) A ( X j(J z ) - J/^(\) + 7 (/i)) , (A2) 

Jz 

where (p(xj) is the propagator of the xj an d A* is the helicity of the photon. In the limit of 
small xj decay width, 

rjsu**'*™ = -"^ (A3) 

In the electric dipole approximation, the xj decay amplitude is written as 

A (xj(Jz) - J/i>W + tM) = E (JJz\L z S z ) A (cc(L z , S z ) - J/^(A) + 7 (//)) 

L Z S Z 



= ^2(JJ z \L z S z )N5 Sz xe* Lz (^ 

L Z S Z 

= N(JJ z \J z -X,X)e Jz , x (fi), (A4) 

where the symbol 5s z \ expresses the heavy quark spin conservation and the normalization 
factor is 

' I ' (A5) 
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Using the results 



Ml — M 2 JU , 
dLips(J/^, 7 )- XJ m 



327r 2 Mj^Ml 



dQ, 



XJ 



'131 



where dQ is a solid angle element in the charmonium rest frame, we obtain 

= Br ( X j - J/^ + 7) E \(JJz\Jz ~ A, A) | 2 ^ / dLips( X j, X)^\A(ij -> xjW + X) | s 

J z 



(A6) 
(A7) 



Br ( XJ - J/^ + 7) E I - A, A) | 2 a(ij - Xj ( J,) + 



This result was also used in the color-singlet-model calculation of Ref . [[7] . 



2. Chromoelectric dipole transitions 



(A8) 



Analogously with eq. (JA5|), we write the cross section of xj production via the chromo- 
electric dipole "decay" of a color-octet ccpSi] state as 

1 -—{o\ar{^ 1 )\i)} 

I 1 nm 



a [ij — > cc 



s sf>\^ X j(J,) + g) = 8(2 j i i )mc {0 ^' {:,Sl ^^ {JJ ^ J ' - s " s * 



X 



M 4 



(A9) 



where (0\Og J ( 3 Si)\0) is a matrix element of nonrelativistic QCD and 



A (ij - cc [ 3 sf \ S z , a])=V2 T a AB Tr 



</*(S z )tf + M) 
2^2 



(A10) 



Oij^lcc is the amplitude for the perturbative process ij — > cc, as given by Feynman rules 



(A, B are the color indices of the heavy quark and antiquark; the heavy-quark and antiquark 
spinors are truncated), and M = 2m c is the mass of the cc state. Combining eqs. ( |A8| ) and 
we have the final result 
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I J — > cc 



J/^(A)+7 + <?) 

Br ( XJ - J/V + 7) Q , OT 1 ^ (O|0 8 XJ ( 3 ^i)|O) E K^l^ " A, *>| S 



xJ2\(JJz\J z -S z ,S, 

s z 



8(2J+ l)m 

2/1 M4 



E |^(u^cc[ 3 5{ 8) ;5 2 ,a])| . (All) 
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FIGURES 

FIG. 1. The Feynman diagrams which describe the leading color-octet mechanisms of J/ip and 
Xj production. The blob represents a nonperturbative transition. The dashed line indicates the 
color-octet intermediate state. 

FIG. 2. Gluon-gluon and quark-antiquark flux factors, <3? 99 and J2q(^gg + ^?g)> plotted as 
a function of r = s/s = x\X% using the GRV-LO parton distributions |l8[ . (a) 7r~-proton reac- 
tions, (b) proton-proton reactions, (c) antiproton-proton reactions. Solid line: <& gg (fiF = M). 
Dotted line: ® gg (nF = M/2). Dashed line: J2q [®qq + $qq] (a*f = M). Dash-dotted line: 
E q [$qq + $qq] (j*F = M/?)- We used M = 3.5 GeV. 



17 



TABLES 

TABLE I. Leading- twist color-singlet and color-octet contributions to J /ip production in ir~N 
collisions at 300 GeV shown together with experimental data |||8). The cross sections have been 
integrated over xp > and normalized by <r e xp 

(all J/tp) = 178 ± 21 nb. The dependence of theo- 
retical contributions on the strong coupling constant a s and on the relative velocity v of the quark 
and antiquark in the bound state is indicated. We used the GRV-LO parton distributions [jig]. 



Observed process 










theoretical subprocesses 


Scaling 


/^r /"oil T /olt\ 


O: 


n~N - 


~> X2 - 


-> J/ip + 7 




U.14o ± U.UzU 






gg - 

qq- 


* X2 

v 3 <j(8) _^ 


a^ 5 


0.059 

0.0041 b 


1.0 a 
0.29 


ir-N - 


-> Xi - 


-> J/V^ + 7 




0.201 ± 0.024 






qq- 


► X15 


a 3 t> 5 


0.0016 


0.19 




qg - 


* X19 




0.0029 


-0.22 




gg - 

qq- 


* Xl5 

► <->i — >• Xl 


a^ 5 


0.0035 
0.0034 b 


-0.27 
0.2 


tt-N - 


-> direct J/ip 




0.56 ±0.03 






gg - 
gg - 
gg - 


♦ ^ - J/V 


a 3 ?; 3 
a 3 t> 7 

3 7 


0.072 
0.16 c 


0.26 
0.5 
1.0 



35(8) 



J/i> 



tt-N ^ all J/tp 



a 3 v 7 



0.020 b 
1 



1.0 

0.028 ± 0.004 



a Reduced to a ~ 0.85 if transverse momentum smearing is taken into account |7j. 
b Using the color-octet matrix elements of Ref. [13]. 



^Combined contribution from 3 -Pj 8 ^ and 1 5'q 8 ^ intermediate states, using the linear combination of 



color-octet matrix elements determined in Ref. [21 
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Fig. 1 



